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Abstract 



We carry out a model-independent EFT method study on the vacuum statistics of general F-term 
SUSY breaking models. Assuming a smooth distribution of Lagrangian parameters, SUSY breaking 
vacua are rare in global SUSY models with a canonical Kahler potential, and have a peaked distribution 
near the cut-off of the SUSY breaking scale in both global SUSY and SUGRA models with a general 
', Kahler potential. After including different mass scales in the Lagrangian, we compare the total number 

(~| ■ of SUSY and non-SUSY vacua and estimate quantitatively the rareness of SUSY breaking. The EFT 

method provides a general view to the amount of parameter tuning needed for a metastable SUSY 
breaking vacuum. The tuning also indicates the importance of R-symmetries in SUSY breaking even for 
metastable SUSY breaking. 



1 Introduction 

Supersymmetry (SUSY) has been proposed for many years to solve several puzzles of the Standard Model. 
To have a physically acceptable model, SUSY needs to be broken dynamically in a hidden sector [3 [51 [31 m [S] . 
Then the mass splitting is transferred to the Standard Model sector by different mediation methods [6l [71 [H [9] . 
The SUSY breaking scale can be anything from a few hundred GeV to the planck scale, depending on models. 
To solve the hierarchy problem and gauge coupling unification, and to provide testable predictions for current 
and near-future experiments, the most interesting models have SUSY breaking at around or not much higher 
than a TeV. 

^ I In the view of the landscape, vast number of models and vacuum states can be constructed from string 

■ " " theory or some other underlying theory [TOl [H [121 [13] , . Each vacuum has its distinct values of physical 

observables, e.g. the SUSY breaking scale, the cosmological constant, etc. If there are enough amount of 
vacua with their observables distributed around the experimental data, one may believe it is very possible 
that at least one vacuum has correct values for all observables, and our universe may live on this vacuum 
through anthropic selection or some evolution processe of the universe. To estimate such possibility, it is 
important to know the vacuum distribution respecting to these observables. One may study a set of models 
build from some underlying theory, e.g. type IIB flux compactification, and hope they have vacuum statistics 
similar to the landscape of the ultimate theory. Many results have been reached in this way [TH [TSl [TB] . 
On the other hand, a simple effective field theory (EFT) with some assumption of parameter distributions 
known from microscopic theories, is useful to demonstrate as good results as more complicated microscopic 
theories can achieve [171 HI] ■ 

We focus on = 4, TV = 1 SUSY models which are favoured by most phenomenology studies. SUSY 
can be broken by F-terms or D-tcrms. The D-term from U(l) vector fields, i.e. the Fayet-Iliopoulos term. 
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is inconsistent with quantum gravity [THl HO] . And non- Abelian D-terms in general can not exist without 
F-terms of a comparable scale, even for metastable vacua [2j. So it is reasonable to concentrate on F-tcrm 
SUSY breaking. The scale of SUSY breaking can be measured by the magnitude of the F-term field strength 



where W is the superpotential for chiral fields {zi}. Note that F has mass dimension 2. For supergravity 
(SUGRA) partial derivatives are replaced by covariant derivatives. 

In section 2, we prescribe the working procedure of the low energy EFT method in a model-independent 
manner. In section 3, taking the superpotential and the Kahler potential to be the most general form, we 
summarize the vacuum statistics in both global SUSY and SUGRA cases. Assuming a smooth distribution 
of Lagrangian parameters, we see that for global SUSY models with a canonical Kahler potential, SUSY 
breaking needs a lot of fine-tuning unless R-symmetries are imposed. In models of global SUSY with a 
general Kahler potential and SUGRA with an either canonical or general Kahler potential, the number of 
SUSY breaking vacua below a cutoff F < Fq is always proportional to Fq, which indicates that the vacuum 
distribution is always peaked at the highest SUSY breaking scale near the cut-off. The distribution respecting 
to the cosmological constant is flat in the SUGRA case. These results are consistent with studies in previous 
literatures [HI [H [21 US] • 

In the string theory landscape studies, one often assumes that there is only one fundamental scale and 
looks for vacua with low energy SUSY. Viewing our set-up as a more general low energy EFT, it is more 
suitable to consider that different parts of the Lagrangian may be generated from different dynamics at 
different energy scales. The presumed mass hierarchy may alter the vacuum statistics. In section 4, we 
scale the superpotential and Kahler potential coefficients to dimensionless by Ms and Mk in addition to the 
Planck mass Mp which enters the SUGRA formula. We find that the vacuum distribution respecting to the 
SUSY breaking scale and the cosmological constant is the same as models with only one fundamental scale. 
The mass hierarchy gives an overall factor to the number of vacua. Comparing to the number of SUSY vacua 
which is obtained using the same EFT method, the overall factor indicates that non-SUSY vacua are rare 
by the factor of Af|/M|-. 

In low energy phenomenology studies, rather than having vacuum distributions. It is more important to 
know the parameter space of models where observables agree with experimental data. The tuning to the 
allowed parameter region could come from anthropic selection, dynamical process at high energy, cosmological 
evolution or just fitting the experimental data. For the SUSY breaking sector, the first step is to estimate the 
size of the parameter space where a metastable SUSY breaking vacuum exists. Our EFT method provides 
a general view to the amount of tuning, which can be obtained from measuring the constrained parameter 
region in the previous statistics procedure. In section 5, we see that in the SUGRA case the superpotential 
gets a large expectation value compared to the SUSY breaking scale, and has to be fine-tuned to get the 
correct cosmological constant. This is a general feature of gravity mediation. A more common tuning of 
order Mg/M^ in both global SUSY and SUGRA cases comes from the metastability condition. Parameter 
region of such a size generally exists for models building. The amount of tuning can be viewed as the allowed 
amount of R-symmetry breaking in models with approximate R-symmetries. This indicates the importance 
of R-symmetries in SUSY breaking even for metastable vacua, as suggested before [MJ [53] . We make final 
comment on the possible extension to multi-field cases. 

2 The general method 

In F-tcrm SUSY breaking models, we have a set of chiral fields {zi\i = 1, . . . ,d}. The superpotential is a 
holomorphic function of {zi}. We also have a Kahler potential K which is a real function of {zi} and {zi}. 
They must be smooth at the vacuum {^^(o)}, so we can expand them near the vacuum in Taylor series. The 




(1) 



2 



most general form is 

W^Y. -J^ r««i-n.(^i - ^1(0))"^ ■■■{zd- MO))"^ , (2) 

^1 • ■ ■ ■ T^d- 

rii 

^= „ 1^ I \ 1^ , Cn,r»i...n,m,(^l - ^1(0))"H^1 " ^l(O))"^ ' ' ' (^d - Zd(0))""(^<i " ^0))"" • (3) 

For general consideration, we include all non-renormalizable terms. Summing indices n's and m's go from 
to cx), and «'s and j's go from 1 to d. To keep K real, the coefficients satisfy 

Cnimi---namd = '^mim ■ ■ -mdnd ' 

We use stars for complex conjugate of coefficients and bars for complex conjugate of other variables. 
For global SUSY, we have the scalar potential 

v^-^i^j.a^^vFa.^w^ (5) 

where K^a is the Kahler metric which is the inverse of the second derivative matrix of K: 

Kjj K^^' = Sf , K^^ = d,^ d,^ K . (6) 
For SUGRA, we have the scalar potential 

1^ = e^(^ KijDs^WD^^W - iWW) (7) 

where the Kahler metric has the same form as ([6|), and partial derivatives in the field strength are replaced 
by SUGRA covariant derivatives: 

F, ^D,^W^ d,^ W + Wd,^ K, D,^ W = {D,^ W)* . (8) 

Metastable vacua are found by searching local minima of the scalar potential V. Right now all mass scales 
in Lagrangian parameters are set to 1, i.e. there is presumably only one scale, and we are searching for vacua 
with F <C 1. We will include different mass scales in later sections. 

2.1 One-field approximation 

The general multi-field model is quite complicated. If there is a field whose mass is much lighter than other 
fields, one can integrate out heavy fields at a low energy scale, or just freeze their values at leading order. 
Then one can have an effective one-field theory which is much easier to study. Fortunately this is the case 
in our consideration. In global SUSY with a canonical Kahler potential, there is always a pseudomodulus 
direction [37]. Its scalar component has zero tree level mass, its fermion component is the massless 
goldstino, and its auxiliary component is the F-term field strength of SUSY breaking. The pseudomodulus 
may get non-zero mass from loop corrections, or equivalently, from non-minimal corrections of the Kahler 
potential [211 ■ The mass is suppressed by the ratio of the SUSY breaking scale to a large mass scale which 
appears in non-minimal corrections of the Kahler potential. The pseudomodulus may also get mass from 
SUGRA interactions. And the mass is suppressed by the ratio of the SUSY breaking scale to the Planck 
mass. Such suppression can be seen in later sections. 

Other than the pseudomodulus, there is no mechanism to introduce a small mass. The natural value of 
field masses should be the same order as of Lagrangian parameters, i.e. of order 1. It is possible by tuning 
parameters to have more than one light fields simultaneously. But such tuning greatly decreases the size of 
the allowed parameter region [17] . For this reason, vacua with multiple light fields only contribute a small 
portion of the landscape unless there is enhancement from symmetries. 
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Now we can write our effective models of the light field z. 

W = ^—an{z - ZqY , K = ^^—Cn^{z-ZQY{z-ZQy, C„m = C„ , (9) 
'-^ n\ '-^ n\m\ 

n n^m 

V=^—dWdW for SUSY, (10) 
ooK 

V = e'^i-^DWDW-3WW) for SUGRA, (11) 
where the d's and D's are respecting to z or z: 

d=—, d=—, DW = dW + WdK, DW = (DW)* . (12) 
oz oz 

2.2 Constraints on parameters 

SUSY-breaking vacua of our interest should have the following properties: 

1. Small SUSY breaking scale, 

2. Stationarity, 

3. Metastability, 

4. Small cosmological constant. 

Each propertie gives some constraint on Lagrangian parameters. In our onc-ficld effective model, coefficients 
an and Cnm can take any values of order 1. But only some region of the parameter space is allowed to have a 
vacuum of our interest. These constraints can be written as a set of Dirac delta hmctions (5's and Heaviside 
step functions 6's. If the vacuum distribution is dense enough on the parameter space so we can take a 
continuous approximation, the number of vacua can be counted as an integral: 

(vacua) = j (i/x(a„,c„TO)2) JJ<5's and 8's . (13) 

The small SUSY breaking scale and small cosmological constant conditions can be expressed as 8(F < Fq) 
and 0(0 < y < Ao) where Fq and Aq are the cut-offs of our interest. For F-term SUSY breaking in global 
SUSY models, the vacuum energy is always related to the field strength as F ~ F^. But the observed 
cosmological constant has a much smaller value. SUGRA coupling or other contributions is needed to get a 
small V. We are not considering those extra effects here. So we just drop off the small cosmological constant 
condition in the global SUSY case. But we include 9(0 <V < Aq) in the SUGRA case. 

The stationarity condition is imposed by requiring the first derivative of V to vanish. This can be 
expressed by the delta function 

5{V') = S{dV)S{dV)\det V"\ (14) 

where V" is the mass matrix 

(15) 



„ _ fddv d'v\ 

{d'^V ddvj ■ 

The form of the delta function gives the correct count of vacua: 

/ <fz6{V') = [ d-zdz ^'(^ - ^(0)) = ^({^(0)}) • (16) 

■' ■' V'(;.(o))=0 

In practice. We treat dV as a function of some coefficient instead of z. The integration of d?ai does not 
give a simple count of aj(o) 's- There is a factor from delta functions: 
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as well as the factor |det V"\ from (|14l) . So we have 

^(n = ^^^^'(a.-«.(0)) . (18) 

The metastability condition is imposed by requiring the mass matrix V" to be positive definite. In the 
one field case here, V" is just a 2 x 2 Hermitian matrix as (ITSt . So the condition V" > is equivalent to 



e{V" > 0) = e{ddV > 0, det V" > 0) . (19) 

With these S and Q functions, we can write down the total number of vacua which satisfies all the 
conditions, 

N{F <Fo)= J d^iSiF < Fo)6{V')e{V" > 0) for SUSY, (20) 

N{F < Fo,0 <V < Ao) = J d^ie{F < Fo)6{V')e{V" > 0)9(0 <V <Ao) for SUGRA. (21) 

As we are to see in the next section, these conditions give constraints on {ao, . . . , a^} if we treat these 
coefficients as variables of S and Q functions. The allowed values of {ao,...,a3} are reduced to small 
intervals. If the distribution of {ao, . . . ,03} is not singular near the allowed region, one can use a uniform 
distribution to approximate the smooth distribution within such small intervals. Thus the integration of 
{ao, . . . , as} picks up a factor proportional to the size of the integration region, which is related to Fq and 
Ao. Other parameters are not constrainec(3. The integration of them just gives a factor which only depends 
on the form of c?/j,(a„, c„m, z). To know the detailed form of d/i one needs to study the microscopic theory 
which produces the landscape. Here we are only interested in the relative distribution respecting to different 
energy scales. Such a overall factor is not included in this work. 



3 Models and vacuum distributions 

We are to apply the method to the general SUSY and SUGRA models. The effective one-field superpotential 
and Kahler potential are given in From now on we set the vacuum at z = by a shift, so we have 

M^ = V^a„z", = V -l-c„„z"z" . (22) 
^-^ n\ ^-^ n\m\ 

n n^m 

At the vcuum, the expectation values for their derivatives are: 

9'WU=o = an, d"d^KU=o = Cnm ■ (23) 

We omit z = in the following notes. Lagrangian quantities are always evaluated at the vacuum. 
In the following models. The F-term i^ 

F = \dW\ = \ai\ . (24) 
So the small SUSY breaking scale condition always gives 



J d^aiQ{F <Fo)^ FdF 



(25) 



We have extracted d^ai from c?/i here since in the small integration region <C 1, the smooth distribution 
is approximated with a uniform distribution. 



^One can also leave an's unconstrained and tune Cnm to fit the vacuum conditions, or spread the constraint in more 
parameters. The difference is just the choice of variables and docs not alter the vacuum distribution result. 

^As we show in later sections, in the SUGRA case, the covariant part of DW vanishes at the vacuum after applying a Kahler 
transformation. So we have F = \ai \ for both SUSY and SUGRA cases. 
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3.1 SUSY with a canonical Kahler potential 

With a canonical Kahler potential K = zz^ the scalar potential is simply 

V = dWdW . (26) 

And its derivatives of our interest are: 

dV = dWd^W = ala2 , (27) 

d^V = BWd^W = alas , (28) 

ddV = d'^Wd^W = ala2 ■ (29) 

Since we are looking for SUSY breaking vacua, we have ai 7^ 0. The stationarity condition requires 
a2 — 0. Then the diagonal element of the mass matrix ddV vanishes. The off-diagonal element d^V also 
needs to vanish to avoid a tachyonic mass eigenvalue. This requires 03 = 0. Now we have a zero mass 
matrix. If non-renormalizable terms are allowed in W , whether the vacuum is metastable is still unknown. 
To insure that there is no high order instability, one can prove that all coefficients in W except ai must 
vanish [23]. Metastable SUSY breaking vacua can only be achieved from a linear superpotential. This is just 
a manifestation of the fact that a SUSY breaking vacuum from any O'Raifeartaigh model with a canonical 
Kahler potential always has a flat pseudomodulus direction j26[ [27] . In our one-field model the only field z 
is the pseudomodulus. The integration intervals for a„'s are infinitesimally small for n > 1, which gives a 
series of infinitesimally small factors to the total number of vacua. For comparison, SUSY preserving vacua 
have ai — which satisfies stationarity and metastability conditions automatically. There is no constraint 
except oi = 0. 

It is well known that SUSY breaking in a global minimum is rare to occur in generic models unless 
one introduces R-symmetries [5U]. Here the result indicates that metastable SUSY breaking is also rare in 
generic models with a canonical Kahler potential. If there is an R-symmetry, usually the pseudomodulus 
z have R-charge 2. To keep W having R-charge 2, all a„'s except oi are suppressed by the R-symmetry, 
rather than by coincidence or fine-tuning in the case without R-symmetries. This suggests the importance 
of R-symmetries for SUSY breaking even in a metastable vacuum. 

3.2 SUSY with a general Kahler potential 

The scalar potential is 



V = ^dWdW . (30) 



And its derivatives of our interest are: 



1 _ _ fr _ _ 

dV = -^—dWd^W - -dWdW 

ddK {ddKf (31) 



= c^la\a2 - c-^^ci2a\ai , 

ddK {ddKY ^ {ddKf {ddKf' (32) 

= Cii a^as - 2c^^ Ci2a^a2 + {Ic^i 0^2 - ^ii cisjaifli , 

2d-dKdd-K (33) 

^ {ddKf {ddKf ' 

= c^la2a2 - c^^cl2a*ia2 - c{^ci2aia2 -t- (2c[f c*2Ci2 - c^iC22)a\ai . 
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Usually the Kiihler potential is written as a canonical form plus corrections, i.e. we set cn = 1 by a field 
redefinition. These quantities can be simplified to 



dV = a\a2 - ci2a[ai , 
d'^V = ajfla - 2ci2ala2 + (2ci2 - Ci3)a*ai 
ddV = a*2a2 — c*2a*a2 — 0120102 + (2c*2Ci2 — £22)0*01 . 



(34) 
(35) 
(36) 



We have oi 7^ 0. The stationarity condition requires 02 = ci2ai. The typical value for c„m is of order 1 
since they are not going to be constrained in our treatment. So both oi and 02 are at the SUSY breaking 
scale oi ~ 02 ~ F. Then we know that the diagonal element of the mass matrix ddV is of order . For any 
values of a„ and c„m, 0202 and 2c|2Ci2a|ai are positive-definite. The rest terms in ddV can be positive or 
negative. So roughly half of the parameter space can make ddV > 0. To make det V" > 0, the off-diagonal 
element d^V also needs to be of order or less, which requires 03 < F. 

Now we have 

03 < oi ~ 02 ~ F . (37) 

Treating dV as a function of 02, we have 



d{dV, dV) 




d{dV) 


a(a*,02) 




da2 



\J\ 

From the estimation of elements of V" , we know 

Idetr' 

So the stationarity condition gives 
d^a2S{V') - 



4a* oi - F^ 



F* 



2 \detV"\ ^2r ^ 

a 02 — (02 - 02(0)) 



1^1 



_p4 

rf^a2-^<5^(o2-02(o))--F2 



And the metastability condition gives 

J d'^as&iV" > 0) - y d^asOidctV" > 0) ^ J 038(03 < F) - 



(38) 



(39) 



(40) 



(41) 



Note the other part of the metastability condition Q{ddV > 0) only reduces roughly half of the integration 
region, thus gives an order 1 factor which we are not interested in here. 
Combining (1^51) (|40p (PT|) we get the number of SUSY- breaking vacua: 



NiF < Fo) 



Ol 



•a3e(F < Fo)S{V')e{V" > 0) 



Fq 



F ■ F^ ■ F^dF - F^ 



(42) 



The power of Fq indicates that the distribution is peaked at the highest SUSY breaking scale. Most states 
live around the cut-off F ~ Fq. In reality, the cut-off is where new physics starts to be important, and the 
low energy EFT becomes not so valid. For example, in the non-SUSY branch of the landscape of type IIB 
fiux compactification [T71 [TS] (which is actually studied as a SUGRA EFT which we are to discuss in the 
next subsection), the cut-off is near the string scale where one has to use the original ten-dimensional theory. 
The vacuum distribution near the cut-off is unknown in the EFT point of view. The distribution from the 
EFT method is more trustable at scales much smaller than the cut-off. 

3.3 SUGRA with a general Kahler potential 

The scalar potential is 

1 



^ " ^ddK^^^^ ~ ^^^^ " e^™(cn («t + a5cSi)(oi + oqCqi) - 3o*oo) . 



(43) 
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These quantities are needed for the convenience of our calculation: 



DW = 


dW + WdK = ai + aoCoi , 


(44) 


dDW = 


d'^W + dWdK + Wd^K = a2 + aicoi + aoCo2 , 


(45) 


d'^DW = 


d^W + d^WdK + 2dWd^K + Wd^K = as + a2Coi + 2aiCo2 + aocos , 


(46) 


Bdw = 


WddK , 


(47) 


d'^DW = 


WB^dK , 


(48) 


BdDW = 


bwBbk + wBb^k . 


(49) 



(51) 



The derivatives of V of our interest are: 

1 - - 3K Bd^K 

= ^"(^^^^^^ + - (»j^)^^^^ - 

= e''""{c:[-^{al + a5c5i)(a2 + aiCoi + aoCo2)+ (50) 
+ (cn coi - c^^ci2){ai + oSc5i)(ai + aoCoi)+ 

- 2oo(oi +O0C01)) , 

1 - - 2BK 2BB^K 

= '''^ddE^^^'^^ ^^BdR- WkT^^^^^^ - WBDW+ 
(BKf _ 2BKBa^K BB'K 2(99^ 

^ BBK BBK {BBKf {BBKY {BBKf ' 
BB^K - 

= e''°°{c^l{al + a5c5i)(a3 + a2Coi + 2aiCo2 + aoCo3)+ 

+ (2c^lcoi - 2c^^ci2)ial + aocSi)(o2 + aiCoi + aoCo2)+ 

- aQ{a2 + aiCoi + aoCo2) + 

+ (cn^cgi + cj"j^^co2 - 2c^^coiCi2 - c^^cu + 2c^^cl2){al + a5c5i)(ai + aoCoi)+ 

- (coi + Cj^i"^ci2)ao(ai + OqCoi)) , 

1 __ _ B^3K - - f)K BFPK 

aif^ii: BKBB'^K BKB^BK B^B^K 2B^BKBB^K -- 
"^^ BBK (BBlCp {BBK^ ^ {BBKf ^ {BBKf ' 
- 2BBKWDW) 

= e""" {c{^ {a*2 + OiCqi + alcQ2){a2 + oiCqi + 00002)+ (52) 
+ (cn^Coi - c{^c\2){a\ + aocSi)(a2 + aiCoi + aoCo2)+ 
+ (cn coi - cj"i^ci2)(ai + aoCoi)(a2 + a^Coi + 03002)+ 

+ (Cn CgiCoi - Cj"]^CoiCi2 - Cj^j^CoiC^a - Cj"]^C22 + 2cj"]^C^2Cl2)x 

X (oi + o3coi)(ai + aoCoi)+ 

- 2ciia5ao) . 

If we apply a Kahler transformation 

W e-^W, K^K + h* + h (53) 
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to set con = c„o = 0, and set the minimal term of K to the canonical form cn = 1 by a field redefinition. V 
and its derivatives can be simplified to 

V ~ a{a\ — Sapao , (54) 

dV — a[a2 — Ci2a\ai — 2a'^ai , (55) 

d'^V = a^aa - 2ci2aja2 + (2c^2 - Ci3)ajai - ala2 - ci2agai , (56) 

BdV — 0302 — Ci2aia2 — 0120102 + (2ci2Ci2 — C22)ajai — 2aQao . (57) 

We have DW = ai at the vacuum after the Kahler transformation. So we can continue using the same 
expression ([25]) for the small SUSY breaking scale condition. The small cosmological constant condition 
requires 

|ai| > V3|ao| > v/|aiP-Ao • (58) 

Since the observed cosmological constant is much smaller than the SUSY breaking scale, i.e. Aq <C F^, we 
have aQ '-^ ai '-^ F. The stationarity condition requires 

02=01201+00^ (59) 
Oi 

which set a2 ^ F. Then we know that the diagonal element of the mass matrix BdV is of order F'^. Roughly 
half of the parameter space can make BdV > 0. To make detl^" > 0, the off-diagonal element d^V also 
needs to be of order F^ or less, which requires 03 < F. 
Now we have 

03 < oo oi 02 F (60) 
Following the same procedure for the SUSY case in last subsection, the stationarity condition gives 

J d^a2S{V') = J d^a2^-^^^SHa2 - 02(0)) ^ J d^a2^d'(a2 - 02(0)) ^ F^ . (61) 

And the metastability condition gives 

J d^asQiV" > 0) - y ^^036(03 <F)^F^ . (62) 

From (|58p we know the integration region of oq is restricted in 

Mei±/F^,.±F).i±iF-^).±F) (63) 

as shown in figure [1] So the small cosmological constant condition gives 

Combining (1^51) (|6ip (p^ we get the number of SUSY-breaking vacua: 

N{F < Fo,0 <V < Ao) ^ J d^ao ■ ■ ■ asSiF < Fo)6{V')e{V" > 0)9(0 < ^ < Aq) 



(64) 



/ F-F^ -F^ ■ AodF - F^Ao 
Jo 



(65) 



Similarly to the result of the global SUSY case, the power of Fq indicates that the distribution is peaked 
at the highest SUSY breaking scale. Most states live near the cut-off F ^ Fq where the distribution is 
unknown in the EFT point of view. The distribution is more trustable at scales much smaller than the 
cut-off. The distribution respecting to the cosmological constant is flat. 
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Figure 1: The value of Oq is restricted in the ring-shaped area in the complex plane to satisfy the small 
cosmological constant condition. The size of the area after including mass scales is also given in parentheses. 



Note that if the Kahler potential is canonical, then 

V = a*ai — SaQOo , (66) 

dV = 2ala2 ~ 2a*Qai , (67) 

d^V = 6a*a3 - 20^03 , (68) 

BdV = 4a^a2 - 2aSao . (69) 

Following the same analysis we still get < oq ~ ai ~ 02 ~ F. The vacuum distribution is the same as 

(ESI). 



4 Including mass scales 

In the previous section we set all mass scales to 1 and look for vacua with <C 1 . This is the case of the 
non-SUSY branch of the type IIB flux compactification landscape where we set the string scale to 1 and look 
for vacua with a small SUSY breaking scale. In phenomenology model building, one usually starts with a low 
energy theory which already have mass hierarchy, and expects some high energy dynamics could generate 
such hierarchy. It would be interesting to see how the statistical analysis of the last section applies to such 
phenomenology-friendly models. Here we assume that there are 3 different mass scales: Ms appearing in 
the superpotential W, which is related to SUSY dynamics; Mk appearing in non-minimal corrections of the 
Kahler potential K, which may, depending on models, come from gauge dynamics, compactified dimensions 
or integrated-out heavy fields; and the Planck mass Mp appearing in the SUGRA Lagrangian. Note that 
W, K, V and z have mass dimensions 3, 2, 4 and 1 respectively. We can write down the general one-field 
SUGRA EFT with these mass scales: 

1^ = V^a„M|-"z", if = V ^c„™M|-"-'"r^'" , (70) 
z — / -J z — / nimi 

71 n,m 

1 - - 3WW WdK 

Now all coefficients a„ and Cnm are dimensionless. We also redefine F to be dimensionless so the SUSY 
breaking field strength is FMg. With the newly introduced scale M5, F is not necessarily small. We are 
more interested in vacua with F of order 1. So the following analysis is a little different than the previous 
section. 
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4.1 Vacuum distributions 



One can do the same calculation as in the last section to get the derivatives of V of our interest. But there 
is a easier way. We take the expressions ([M)) -([57 1) . do the following replacement: 

(72) 

and insert powers of Mp to where the mass dimension does not match. Then we get 

A,f 6 

V = alaiMl - 3a*aoj^ , (73) 



r2 



dV = ala2A4 - ci2at ai ^ - 2a* a, ^ , (74) 



ddV = ala2Ml - (0^20^02 + 0120102)—^ + (2ct2Ci2 - 022)0^01— f - 20000-7! . (76) 

lVl}i -'"AT iVlp 



There is one issue we would like to clarify before we continue. The results (|54l) - ([57|) and ([75 )) - ([75)) have 
been simplified by several field redefinitions. In the previous section all mass scales are set to 1. After field 
redefinitions, values of Lagrangian parameters change amount of order 1 and no singularity is introduced to 
the parameter distribution. So one can continue the analysis using the newly defined fields. This is not so 
obvious with mass scales included. For the translation (z — zo) — > z, the typical value of zo is of order M5. 
After the translation, values of a„'s change amount of order 1, and values of c„m's change amount of order 
MsIMk- The Kahler transformation to set co,i ~ Cno = is 

W^e^^W, K^K + h* + h, /i = - V i-co„Af2-«2" (77) 

n 

where we have included mass scales and the required form of h is given. Giving the condition Mx < Mp, 
values of a„'s change amount of order M^/Mp at the leading order. Finally, the field redefinition to 
set cii — 1 scales all coefficients by factors of order 1. In summary, values of dimensionless Lagrangian 
parameters change amount up to order 1 after all field redefinitions. If we start with a smooth distribution 
of Lagrangian parameters, the distribution after field redefinitions is also smooth and can be approximated 
by a uniform distribution. So it is valid to use simplified results ([75|) - ((75)) for distribution analysis. 

We still have F = \ai\ after the Kahler transformation. So the small SUSY breaking scale condition gives 
the same factor as ([25]). The small cosmological constant condition requires 



KI>^|a„l>^KP-^. (78) 

Since the observed cosmological constant is much smaller than the SUSY breaking scale, i.e. Ao ^ F^M|, 
we have 

Mp , , 

ao-F— . (79) 



The stationarity condition requires 



Ms , , ,ai Ml , 1 , 1 



a2 = e.2ar^ + ^a^^^ ^ FMs{— + — ) . (80) 
Then we can estimate the magnitude of the diagonal element of the mass matrix: 

aay^F^M^(^ + ^). (81) 



11 



From now on, conditions Ms ^ Mk and Ms ^ Mp are imposed to keep only the lowest order terms in 
estimations. Although in most realistic models Mk is either much smaller than or at the same order as 
Mp, their magnitude relation is not fixed right now (and is to be studied in later subsections). So we keep 
both lowest order terms in the result. Roughly half of the parameter space can make ddV > 0. To make 
det V" > 0, the off-diagonal element d^V also needs to be of the same order or less. We have 



d^V ^ a.FM^s + F'Mli-^ + TTt) ^ ^^^^siTIT + TTt) > (82) 



'Ml My- ^^M^ Ml 



a3<FM|(^ + ^). (83) 



which requires 



The field z has mass dimension 1, which makes the previous count of vacua (1161) having mass dimension 
-2. To make a dimensionless count one needs to rewrite the delta function as 

6Hz - Z(o)) ^ '52(^(z - Z(o))) = Ml5^(z - z^o)) ■ (84) 



So S{V') picks up a factor of M| compared to ([14]) and (fT8| : 

SiV) = MlS{dV)5idV)\detV"\ = Ml^-^^^5\a, - a,(o)) . (85) 
Treating dV as a function of 02, we have 



\J\ 



didV, dV) 



From estimations of elements of V", we have 



d{dV) 



da2 



2 



AalaiM^ - F'^M^ . (86) 



l<ietrVF.M|(-^ + -^), (87) 



So the stationarity condition gives 

J £a,5{V') = J d^^Ml^-^^S'ia, - a,^o)) - F'^U^ + ^) . (88) 
And the metastability condition gives 

/ d^aMV" > 0) ^ / d^aMa. < FMUj^ + ^)) - F^M|(^ + ^) . (89) 
From (j78p we know the integration region of aq is restricted in 



Mp / Ao Mp^. , Mp Ao.Mp 

which is a similar ring-shaped area as shown in figure [TJ So the small cosmological constant condition gives 

9 . , ^ FMp AqMp AqMI , , 

<i%e(0<V<A.,^2,^.^^~^^ (91) 
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Combining (gH]) ([SS]) (IMl) (EB) we get the number of SUSY-breaking vacua for the general SUGRA case: 



N{F < Fo,0 < V < Ao) / d^ao ■ ■ ■ aaQiF < Fo)S{V')e{V" > 0)6(0 <V <Ao) 



F.F^MU^ + ^)-F^Mli^ + ^)--^^dF (92) 



~ F^AnMlMU-^ + -^) 
s P'^Ml^ Ml' 

In most realistic models Mk ^ Mp, so we have 



N{F <Fo,0<V<Ao)^ Po^^^^^f^ ■ (93) 



Comparing to the result ([65)) before including mass scales, the vacuum distribution has the same form 
respecting to F and A. The power of Fq indicates that the distribution is peaked at the highest SUSY 
breaking scale. Most states live around the cut-off F ^ Fq. The distribution respecting to the cosmological 
constant is flat. But now we have the extra factor from including mass scales. With the mass hierarchy 
presumed, we have made our analysis procedure valid up to i^o 1- The distribution is reliable all the way 
up to the cut-off. But the tail of the distribution after the cutoff could still be complicated. 



4.2 The Global SUSY and canonical Kahler potential limit 

If we take the limit Mp — > oo of ([5^ . the number of vacua goes to oo. This is because there is a factor coming 
from the small cosmological constant condition which set the value of oq having an allowed area proportional 
to Mp. As we have discussed before, the small cosmological constant condition should be dropped off in 
the global SUSY case. Most analysis goes through similarly to the SUGRA case. Combining (|25p([88|)((89|) 
and taking the limit Mp — >■ oo, we get the number of SUSY-breaking vacua for the global SUSY case with 
a general Kahler potential: 

d\r ■ ■ ■ a,e{F < Fo)5{V')Q{V" > 0) - F- F^J^ ■ F'j^dF ^ F^j^ . (94) 

Comparing to the result (j42p before including mass scales, the vacuum distribution has the same form 
respecting to F. But now we have the extra factor from including mass scales. 

To see the vacuum distribution with a canonical Kahler potential, we take the limit Mx — oo, (|92p 
becomes 

A 

N{F < Fo,0 <V < Ao) ^ F^^^ . (95) 

The distribution respecting to F and A has the same form as ([92]) . This explains what we have seen in the 
SUGRA case before including mass scales where we get the same distribution with an cither canonical or 
general Kahler potential. 

If we take the limit AIk — oo of ([M)) . the number of vacua goes to 0. This shows that metastable SUSY 
breaking is rare in generic global SUSY models with a canonical Kahler potential unless one introduces 
R-symmetries, as we have discussed before. 



4.3 Comparing to SUSY vacuum distributions 

Since we only assumed a uniform parameter distribution in general cases, the vacuum distributions what 
we have got before have only relative meaning. Even though without knowing the the detailed form of 
dfi{a„, Cnrm z) ffom the microscopic landscape, the factors we have got from including mass scales enable us 
to do a general comparison between the distributions of non-SUSY and SUSY vacua, which is an interesting 
topic to study in the landscape. 
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We start again from ([751) - ([7^ . For SUSY vacua, the F-term vanishes, so ai = 0. Then we always have 
< 0. In the SUGRA case it is still worth to count the vacua with — Aq < < and expect that other 

corrections of the same order could set the cosmological constant to a positive value. To get such a vacuum 

requires 

dV always vanish for ai — 0, so the stationarity condition gives no extra constraint. The elements of the 
mass matrix become 



d V = -a,a,j^^ ^ a,^^ , (97) 



BdV ^ a^azAfJ - 2aQaoj^ - ala2M^ . (98) 

The dominant term of ddV has a positive definite form and \d'^V\ ^ ddV. detV" > is always sati- 
fied except in an exceptionally small parameter space region where 02 is of order \/T^/ {MsMp). So the 
metastability condition also gives no extra constraint. 

There is no constraint on a„ for n > 1. They can take any value up to order 1. It is ai = that satisfies 
the stationarity condition. So we treat dV as a function of ai and a^. Then we have 



1^1 



d{dV,dV) 



- a*a2Af| - 4a5ao-^ ^ M| . (99) 

Ivl p 



d{a\, ai 

From the estimation of elements of V" ^ we have 

A 

|det V"\ ^ {a^a^fMl - «2«2^^ ^ Ml . (100) 

So the stationarity condition gives 

J d'a.SiV, F = 0) = J d\,M's^-^^S'{a,) ^ 1 . (101) 

The metastability condition gives only an order 1 factor. From (|96l) we know the small cosmological constant 
condition gives 

/A Af'^ 
d^aoe{-Ao <V<Q)^ . (102) 

Combining ([TOT|) (fTOS ]) we get the number of SUSY vacua for the SUGRA case: 

/A Af^ 
d^ao ■ ■ ■ a^SiV, F = 0)Q{V" > 0)6(0 <V<Ao)^ -j^ . (103) 

For the global SUSY case, most analysis is similar to the SUGRA case. Dropping off the small cosmological 
constant condition and taking the limit Mp — >■ 00, we get the number of SUSY vacua: 

N{F = 0) - / d'^ai ■ ■ ■ a3(5(U', F = 0)&{V" > 0) - 1 . (104) 



Now we can compare these results to (|93l) and ([94]). Setting Fq ^ 1, wc have 



for SUSY (105) 
mF^l,0<V<Ao) ^^^^^^^^ 

7V(F = 0, -Ao < U < 0) M|. ^ ' 

In both SUSY and SUGRA cases, if all mass scales are set to be 1, the numbers of SUSY and non-SUSY 
vacua are comparable when the cutoff Fq is set to be of order 1. If we assume Ms ^ Mk for low energy 
model building, non-SUSY vacua are always rare compared to SUSY vacua by the factor of Af|/M|.. 
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5 Parameter tuning 



In previous sections we have stated many times that the result indicates the rareness of SUSY breaking. At 
the end of the last section we even get the rareness quantitatively. These results come from the assumption 
that Lagrangian parameters have a smooth distribution and can be approximated by a uniform distribution. 
Such an assumption can not go beyond the tree level SUSY breaking branch of the landscape. In the 
intermediate or low energy SUSY breaking branch where SUSY is broken dynamically, there are singularities 
in the parameter distribution, and the vacuum distribution regarding to F is very different than the tree level 
SUSY breaking branch [30l [311 [32] . Besides the counting of metastable vacua, non-perturbative instability 
may furthermore alter the distribution [33 , 34 , J51 [33] ■ Cosmological settings such as chaotic inflation also 
realize vacua with different probabilities ^37. .38. ,3^ . For low energy model building, it is more important to 
know which parameter region leads to phenomenologically acceptable models rather than having the vacuum 
distribution. One can ask how much tuning is needed to get a desired vacuum with appropriate F and A, 
and leave the question of the origin of the tuning to the future study of fundamental theories. Although it 
is quite a trivial routine to work out the allowed parameter region from a specific model, our EFT method 
provides a general view to the amount of tuning. In previous sections, constraining Lagrangian parameters is 
a procedure to calculate the vacuum distribution. But the constrained region of a„'s is just what we concern 
here. 

We assume Ms ^ Mk < Mp and F ~ 1. In the SUGRA case, (|90)) shows ao is restricted in a ring-shaped 
area as shown in figure [1] with its magnitude around 

Mp Mp 
V3M5 Ms 

The large value of ao seems unnatural, but is required to cancel the contribution from DW if there is no 
other contribution to the vacuum energy. SUSY breaking gives the gravitino mass 



m3/2 = ^-ao^~-^ (108) 



W _ Ml Ml 

mI^^^^mI ^ M~p 

which means that SUSY breaking occurs at an intermediate scale between 7713/2 and Mp. This is actually a 
feature of phenomenological SUGRA models, i.e. gravity mediation. As a comparison, other scenarios such 
as gauge mediation favor low scale SUSY breaking. 

In the previous treatment, a„'s are scaled to dimensionless by powers of M5. We have assumed that 
typical values of a„'s are of order 1 if they are generated by the same dynamics which breaks SUSY. But 
ao does not contribute to SUSY breaking except for adjusting the cosmological constanl[f|. It is reasonable 
to think that oq gets contribution from processes other than the SUSY breaking dynamics and not related 
to Ms- Alternatively, one can introduce another hidden sector which dynamically generates ao with its 
magnitude to the correct order, as what has been done in retrofitting models [40 . However, the narrow 
thickness of the ring-shaped area as shown in figure [1] suggests that an order Ao/A/| tuning is necessary for 
a small cosmological constant. 

Now we focus on constraints on other parameters which are common in both global SUSY and SUGRA 
cases. We have seen |ai| = F which is just the SUSY breaking field strength. The value of a2 is determined 
from oi by a delta function. This is just the requirement of a stationary vacuum and should not be viewed 
as tuning. In fact, when we do the translation (z — zq) — z to set the vacuum at z = 0, the value of zq 
occupies two dimensions of the parameter space which compensate the lost two dimensions from fixing a2. 
The only tuning comes from the metastability condition. The constraint (|83p implies 

Since takes a complex value, the needed amount of tuning for metastable SUSY breaking is of order 
Mg/Mj^. Conversely, the metastability condition for SUSY vacua is always satisfied in the global SUSY 



^Even ao appears in DW in the SUGRA formula, it is suppressed by powers of Mg/Mp or Mjf/Mp. 
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case and only excludes exceptionally small parameter region in the SUGRA case. No tuning is needed to 
get a SUSY vacuum. Note that we have done several field redefinitions to simplify the form of the potential. 
The relation between a„'s and parameters of the original model could be complicated. So the amount of 
tuning may be spreaded in several of the original parameters when we try to build a realistic model. 

Our analysis up to now is based on tree level calculation. One may argue that loop corrections could 
have extra contribution to the mass of z, thus relaxing the tuning of (I109p . One example comes from direct 
mediation where the SUSY breaking field z couples to some messenger fields 4>i in the same hidden sector. 
The pseudomodulus gets a Coleman- Weinberg mass of order M'^/M^ where is the typical messenger 
mass . The Coleman- Weinberg formula is only valid when ^ Ms- So the contribution to the mass 
of z is much smaller than Ms- One still needs some tuning for 03 to keep the vacuum metastable. In fact, 
one can integrate out heavy fields when Ms and get non-minimal corrections to the Kahler potential. 

Then the Mk can be substituted by M^ in the tuning of (|109p . 

Although much tuning is needed when Ms <^ Mk, there is always some allowed parameter region as 
long as non-minimal corrections of the Kahler potential do not vanish. Precisely speaking, if we follow the 
procedure how we get (|83p . the leading order contribution to (|109p is from c^2 and C22- If one can freely 
choose the form of the superpotential and the Kahler potential, non-vanishing C12 or C22 generally allows 
one to build models with metastable SUSY breaking vacua in some parameter region. This feature has been 
observed previously while building models locally equivalent to the Polonyi model |42) . Here our result shows 
that the argument is true for most general F-term SUSY breaking models in both global SUSY and SUGRA 
cases. The size of the allowed parameter region is always proportional to Mg/Mf.. 

We have not introduced any symmetry in our analysis. In SUSY breaking model building, R-symmetries 
are often applied because of their connection to SUSY breaking global vacua dS] ■ Usually the pseudomodulus 
z has R-Charge 2. z can only appear linearly in the superpotential because of the R-symmetrjQ. Although 
only part of the parameter space can make a metastable SUSY breaking vacuum, the reduction is usually 
of order 1 and no tuning is required. Once the R-symmetry is broken, terms like a^z^ may appear in the 
superpotential and destroy the metastability. So the tuning of ()109|) also characterizes the allowed amount 
of R-symmetry breaking which can be introduced into R-symmetric models. We see again the importance 
of approximate R-symmetries in metastable SUSY breaking model building, as discussed in literatures for 
different reasons [24l [25] . 

Finally, we would like to comment on extending the one-field approximation to multi-field cases. We only 
studied the behaviour of the potential along the pseudomodulus direction in our current analysis. Other 
fields, if they are heavy than the SUSY breaking scale, could be viewed as appearing in the non-minimal 
corrections of the Kahler potential after being integrated out. If more than one fields are light, one can 
still identify the pseudomodulus-goldstino along the F-term direction. The same tuning of (|109p is needed 
following the same analysis of the one-field case. In addition, stabilizing other field directions may further 
reduce the allowed parameter region, either by a factor of order 1 or more tuning. Having more than one light 
field seems statistically rare, but may be preferred by symmetry consideration or for solving phenomenology 
problems 03111113]. 
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